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ANALYTICAL TRIPPING LOADS FOR 
STIFFENED PLATES ’ 



D. A. Danielson 
Mathematics Department 
Naval Postgraduate School 
Monterey, CA 93943 

January, 1994 



ABSTRACT 

The subject of this paper is the buckling behavior of a rectangular plate, with parallel thin- 
walled stiffeners attached to one side, subjected to a combination of axial compression, lateral 
pressure, and bending moment. The plate is modeled by the Von Kdrmdn plate equations and 
the stiffeners by a nonlinear beam theory recently derived. An analytical solution is obtained 
for the buckling load corresponding to a torsional tripping mode of the stiffeners. The effects 
of various boundary conditions, imperfections, and residual stress are included. 

1. INTRODUCTION 

Stiffened plates are a basic structural component of ships and submarines. These struc- 
tures are designed with generous safety margins against overall collapse triggered by buckling. 
The object of analytical work is to determine design criteria to inhibit buckling at any stress 
less than yield. Recently [see Danielson et a/ (1993)], we have developed an analytical formula 
for the buckling load of a stiffened plate subjected to a combination of axial compression 
and lateral pressure. The object of the present paper is to improve and extend our previous 
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analysis. A review of the literature, given in our earlier work, will not be repeated here but 
details of the analysis, which supersedes our earlier work, will be recorded here. 

We first consider a plate which is initially rectangular in shape and has several parallel 
I-stiffeners spaced a distance b apart. The structure is subjected to a combination of uniform 
axial compressive stress a (force per unit area of a side), uniform lateral pressure p (force per 
unit lateral area of the plate), and uniform bending moment M (moment per unit length of 
an edge). We suppose that at low values of <7, p, and M the plate and stiffeners simply bend 
and compress symmetrically. Our object is to find the critical load at which the stiffened 
plate may buckle into an alternate mode (see Fig. 1). 

Our present analysis is based on the following simplifying assumptions; 

(i) Each plate-stiffener unit of width 6 undergoes an identical deformation. 

(ii) The plate obeys the nonlinear Von Karman plate equations [see Timoshenko and Gere 

(1961)]. The stiffeners obey the nonlinear beam equations derived by Danielson and 
Hodges (1988). 

(iii) The plate and stiffener material is elastic, linear, and isotropic. 

(iv) Every particle on the bottom surface of a beam undergoes the same displacement as 
the corresponding particle on the top surface of the plate, and every line of particles 
in the beams normal to the plate surface remains normal to the deformed plate at its 
surface. In other words, the bases of the stiffeners are clamped to the plate. 

(v) The prebuckling displacements are less than the maximum thickness of the structure 
and independent of the transverse coordinate. 

(vi) The incremental buckling extensional strains at the midsurface of the plate are negli- 
gible. 
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(vii) The incremental buckling displacements may be approximated by the fundamental 
harmonic in their Fourier expansions. 

(viii) The plate and beams are so thin that their thicknesses are negligible compared to 
their width, height, length, and the wavelength of deformation. A stiffener is so slender 
that its width and height are negligible compared to its length and the wavelength of 
deformation. 



2. GENERAL POTENTIAL ENERGY FUNCTIONAL 

It follows from assumption (i) that we need only analyze a single plate unit containing 
a single stiffener. From assumptions (ii)-(iii), the potential energy of the plate plus beam is 
given by: 

P[u,u,tn] = /I ^ + Y + (1 - i^)en] 

+ — pw -f Mwii > dxi dx2 



+ D ^ + UWUW 22 + ^ + (1 - l')wl 



I J L 26*712 + 26*713 j dx2dx3 + aA{ui - cwu)dxi 

(1) 

Here (xi, x^, Xj) are Cartesian coordinates measured from the midpoint of a side of the plate. 
The plate unit has length a, width b, and thickness /, while the beam has cross-sectional 
area A and centroidal height c. The elastic constants are defined by 



D = 



Et^ 



G = 



12 ( 1 - 1 / 2 )’ - 2{\Av) 

where E is the Young modulus and t/ is the Poisson ratio. The displacements of the plate 
midsurface in the xi^X2^ and 2-3 directions are denoted by u{x^, X2), v(x],X2), and 2'2), 

respectively. Subscripts on u,n, or w denote partial differentiation with respect to the 
coordinates Xi or 2:2; e.g., W12 = The extensional strains at the midsurface of the 

plate are given by 
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en = ui + -tv^, ei2 = -(u2 + + W1W2), 622 = ^2 + 2^2 

The strains in the beam are denoted by 711(11, X 2 , 2^3), 712(^1? ^21 ^ 3)1 713(11, 2^21 ^3) 

and are related to the displacements by equations (9, 10) of Danielson and Hodges (1988), 
which upon invoking eissumption (iv) are transformed into: 



7 11 = ■£'11 + (■£'12 + Ci 2)<^3 ~ ■£' 13*^2 + 2^2 "h 2*^3 

712 = E \2 — -■£' 11 <^ 3 , 713 = ■£'13 + 

■£^11 ~ ^11 ~ ^ 3^11 + -^^ 112 i ■£'12 ~ ^12 “ ^^ 3^12 + 2 '^ 2^12 

^ 1 _ 1 , _ 1 , 

E\z — -^21^12 + 2 ^^ 3^12 + 2^W\\W\2 

, 1 _ 1 , _ 1 , 

92 — ~ 2 ^ 2^12 + 2 ^^ 3^12 ~ 2 ^‘’^ n ’’^^2 

■ 4 >z = “ 2 ^ 3^12 



Here A(i2, 13) is the Saint-Venant warping function for the beam cross section; subscripts 
on A denote partial differentiation with respect to X 2 and 13. Bars over a symbol denote its 
value at the beam axis; e.g., the axial displacement of the beam centroid is u(ii) — cuJi(ari). 

Substituting these relations into (1) and neglecting higher than cubic terms in the dis- 
placements (these are not needed in our subsequent analysis), we obtain a lengthy expression 
for the potential energy which forms the basis for our subsequent analysis. Among all the 
functions satisfying the geometric or natural boundary conditions the one which causes the 
potential energy to be a minimum is the equilibrium state. We suppose that the outer edges 
of the plate are free to displace in the horizontal plane but restrained in the vertical direction, 
so the geometric boundary conditions are 



u;(0,a'2) = w{a,X2) = 0 

Note that the case of simply supported edges is obtained by setting M = 0. 



(2) 
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3. PREBUCKLING SOLUTION 



The prebuckling equilibrium state is denoted by {u,w). It follows from assumptions 
(ii)-(v) that the potential energy in the prebuckling state is 



Jo 



E{A + ib) Sj {bD + E/22) ° A 

1 r^/ic 

+a(A + tb) Ui —aAc u>u —pb w -\-Mb Wnj dxi 

Here I 22 is the moment of inertia of the beam section about the 



(3) 



hi = / L.n, xldx2dx3. 

J J ttzxion 



The prebuckling displacements are determined by the variational equation 



6P = 0 



Taking the variation of (3), and integrating by parts with respect to ij, we obtain the 
differential equations 

(6/) + E/ 22 ) u?iiii — Ei4c uiii — = 0 . 

(4) 

— EAc ^11 +E(j4 + tb) Uj Acr{A A tb) = 0 

and boundary conditions 

at X = 0 and x = a : {bD + E/ 22 ) ^11 —EAc Ui —aAc A Mb = 0 (5) 



The solution to the linear boundary value problem (2), (4), (5) is 






0 

W 



a Ac Wii 
E AAtb 



pxija - xi)(g^ + axi - x^) 

p 42-2 

24(E + -(/22- j^)] 



Mx\{a — xi) 



2[E + ~r{hi — 






A "I" tb 



)] 



(6) 

(7) 



Note that (6)-(7) reduce in the case /22 = = 0 to the well-known exact solution for an 

isolated wide plate, and in the case t = 0 to the well-known exact solution for an isolated 
beam. 
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4. BUCKLING SOLUTION 



According to the energy criterion of elastic stability, the prebuckling equilibrium state is 
stable if and only if the energy functional which represents the increase of the total potential 
energy in a displacement field to some slightly adjacent state (u; +u;) is non-negative: 

P\w -ftxj] — > 0 . (8) 

Since the prebuckling state is an equilibrium state, the terms in (8) which are linear in the 
incremental displacement w must vanish. It follows that the terms Q[u^] in (8) which are 
quadratic in the incremental displacement must be non-negative; 

QM > 0 . 

The critical case of neutral equilibrium occurs when there exists a buckling mode w„ satis- 
fying 

(^[u^cr] = 0 (9) 

Q[w 7 ^ tCcr] > 0. (10) 



The eigenvalues <Tcr,Pcr, and Mcr which render (9) zero are the critical buckling loads. 

From the first integral in (1) and assumption (vi), the quadratic functional for the plate 



IS 



»— 4f{ 



D 



^11 ^22 / , ^ 
— ^ h VW\\W22 + 1" (1 ~ 






dX2 . 



(11) 



From the remaining integrals in (1), the quadratic functional for the beam is 
EHi 






El Ac 

/ w\^2dxi-)f / tx;]] -1-^ / w\^dxi 

Jo 2* Jo 2 Jo 

(12) 



Here 1 is the polar moment of inertia about the Xj-axis, J is the Saint- Venant torsion 
constant, and /fj, H 21 Hz are constants defined by the following integrals over the beam 
cross section: 
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I 

J 

H, 

H2 

H, 



I /b.*m xl)dxidx3 

J •/•ectioD 

/ /b..m 1(^2 + -^ 3 )^ + (X3 — ^2)^]dx2 dx3 

J •/•eclion 

/ /b..m \'^dx2dx3 

J */«ection 

/ /b.»m ^ 3 ( 3^2 + 3 :^ - ^2 - ^ 3 ) + 2A(l2 + A3)]<^2-2 dx3 

J •'•ection 

7 / /be»m ^3[(3T2 “ -^ 3 )^ + (^3 + -^ 2 )^ + + ^3 ~ -^2 ~ ^3)]^^2 ^^3 

4 •/ •/•ection 



The total quadratic functional for the plate plus the beam is the sum of ( 10 ) and ( 11 ). 

Next, we calculate the cross section properties for a beam composed of a thin web and a 
thin bottom and top flange. The web has thickness and height h^; the bottom flange has 
thickness tb and width /ij,; the top flange has thickness tj and width hj. The Saint- Venant 
warping function for this thin-walled cross section is: 



A = ^ 



X2{2h^ + tj - X3) flange 

122:3 web 

—X2X3 bottom flange 



Using approximation (viii) we obtain: 



A 


= 


ib^b ^ i XI.' h yj ^ i J h J 


Ac 


= 




I22 


= 


t 


I 


= 


h^b tjhj , . 1 2 1 

12 + 3 +12 


J 


= 


tih^ t)h, 
3 3 3 


H, 


— 


ifhih) 

12 


H2 


— 


t)Kh, tlhl 
3 6 


Hz 


= 


1 / ;,3 J, 1 
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In accordance with approximation (vii) and the boundary conditions (2) and (5), we 
represent the incremental buckling displacement by the following shape (an arbitrary multi- 
plicative constant has been set equal to 1): 

(13) 



. mirxi . irx2 

ly = sin sm — , m = l,2,3 

a 0 



Substitution of this buckling mode into (11) plus (12) and application of the inequality (10) 
leads finally to: 






A tb 



J A 

- -T-r-r.W. 



F 42-2 



_ A^<? . 



hD ,mb 
T a 



t>JJ ,m0 a 

+ — ( — + — r) + 

mo 






T 



(14) 



Here m is taken to be the integer which gives the lowest value of O’er in (14). Note that 
(14) reduces in the case tj, = = ty = 0 to the well-known exact solution for an isolated 

plate, in the case per = Mcr = t = 0 to our previous exact solution for an isolated beam [see 
formula (28) of Danielson et al (1990)], and in the case Per = = 0 to formula (87) (with 

^'pbe — (*^e)^e ^nd I -f P Rs H\) of Adamchak (1979). 



5. CLAMPED EDGES 



In this section we consider the clamped case when the bending moment M is not pre- 
scribed, but the rotation at the edges is completely restrained. Then the additional geometric 
boundary conditions are 

iui(0,i2) = ( 0 , 2 : 2 ) = 0 (15) 



Note that in this case each material particle on an edge cross-section is totally restrained 
from any motion at the buckling point. 

The prebuckling solution to the linear boundary value problem (2), (4), (15) is (6) and 

pxl(a - ii')2 



O 

W= 



F 



(16) 



-f Gj 
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Note that (6) and (16) reduce in the Ccise I 27 = .4 = 0 to the well-known exact solution for 
an isolated wide plate, and in the case t = 0 to the well-know'n exact solution for an isolated 
beam. 

In accordance with approximation (vii) and the boundary conditions (2) and (15), we 
now represent the incremental buckling displacement by the following shape; 



m-KXx . TtX2 ... 

lu = (1 — cos )sin— m = 2,4,6, 

a b 



(U) 



Substitution of this buckling mode into (11) plus (12) and application of the inequality (10) 
leads finally to 






<^cr < 






bD.mb a , 2a^ 

+ ^l(— + —J + — 



' 2 a ' mb' ' m'^b^ 



+ GJ 



/ + 



2;r2 



(18) 



Note that (18) reduces in the case pcr = t = 0 to our previous exact solution for an isolated 
beam [see formula (28) of Danielson ei al (1990)]. 



6. IMPERFECTIONS 

In this section we suppose that the structure has an initial normal deflection in the shape 
of the prebuckling normal displacement caused by the pressure p. Specifically, w'hen M is 
prescribed on the edges the initial normal deflection is 

16 ari(a - Xi)(a^ + axi - x\) 

5^5 

while for clamped edges the initial normal deflection is 

16 j?(fl - 
a* 

We also assume that the amplitude W of the initial displacement is less than the maxi- 
mum thickness of the structure. Then the prebuckling displacements are still given by our 
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previous solutions (6)-(7) or (16), and the quadratic functional Qpiate is still given by (11)- 
The only effect of this prebuckling deflection is to create a new term in Qbeam which is the 
same as the middle integral in (12) with the prebuckling normal displacement replaced by 
the initial displacement, and vfe can use our previous calculations to evaluate this integral. 

We thereby find that the effect of this imperfection is to add an additional term to our 
previous formula for (Tct' 

M case: 



t I Ar 






A tb 



Ccr = cr„[w = 0] + 



5q2 



, th^ 



Clamped case: 



(19) 



Oct = Oct [IV = 0] + 



I 






( 20 ) 



Note that for an asymmetrical structure torsional deformation of the stiffeners may ini- 
tiate upon application of the slightest load, so bifurcation may not be able to be used as the 
buckling criteria [see Ostapenko and Yoo (1988)]. 



7. RESIDUAL STRESS 

The simplest way to account for residual stress is to assume that the plate is subject to 
a uniform compressive residual stress 5, while the beam is subject to a counterbalancing 
distribution of residual stress <7r(i2>^3) [see Hughes (1983)]. The only effect of this residual 
stress is to create a new term in (?pu„ which is the same as the last integral in (11) with 
<7 replaced by S. Note that the analogous term in Qbeam is zero because iZJ(ari ) = 0 for the 
assumed mode shapes (13) and (17): 

-III ^dx,dx,dx,=^0 

J J «/be«m ^ 
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We thereby find that the effect of this residual stress is to add an additional term to our 
formulas for Ccr- 

<7cr = = 0] (21) 

^ ^ ~w 

CONCLUSIONS 

Simple analytical formulas that include the effects of combined loading, various boundary 

conditions, imperfections, and residual stress do not appear to exist in the literature. As 

an example of the numerical predictions of our formulas, let us assume the following typical 

parameter values [taken from Smith (1975)]: 

E — 30,000 ksi <6 = 0 
i/ = .3 /?6 = 0 

a = 48 in <u, = .28 in 

6 = 24 in h^, = 5.5 in 

< = .31 in </ = .56 in 

hj — 3.1 in 

For this example, the formulas (14), (19), (21) and (18), (20), (21) reduce to 
M case (m = 1): 

= 49 + 60p„ + AbM„ + 55W - .765 (22) 

Clamped case (m = 2): 

= 105 + 6.6per + 30VF-.765 (23) 

The collapse loads of <Tct = 27.8 psi {pcr = 0) and <Tct = 27.1 psi {pa- = -015 ksi) measured 
by Smith (1975) on ship grillages may be accounted for by choosing appropriate values of 
Mcr,W, and 5 in (22) or (23). 

We have made an attempt to verify the accuracy of some of the approximations upon 
which our analysis is based. For instance, in the M case, we included the effect of beam 
cross-sectional deformation by allowing the web to undergo a lateral buckling displacement 

V = sm (C1X3 -f C2X3) , 

a 
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where C\ and C 2 are constants determined by minimizing the critical axial stress <Jct- For 

the parameter values listed above, this effect turned out to be negligible. It is possible to 

invent unusual cases in which the cross-sectional deformation is of importance, but for most 

» 

practical dimensions the assumption of a rigid cross-section seems ok. 

For another instance, in the clamped case, we added together the m = 2 and m = 4 

normal buckling displacements (17): 

, 27TI1 47rJi . 7TX2 

ty = [1 — cos 1- 63(1 — cos )J sin — r— , 

a a b 

where C3 is a constant determined by minimizing the critical axial stress Cct- For the 
parameter values listed above, this effect also turned out to be negligible. For most practical 
dimensions the assumption of a simple buckling mode seems ok. 

At any rate, if we include the additional displacement functions needed for the above 
effects, it doesn’t seem possible to obtain simple formulas for the unknown coefficients or 
buckling loads. 
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